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Abstract 

We consider a generalisation of the vicious walker problem in which A'^ random walkers in R'* 
are grouped into p families. Using held-theoretic renormalisation group methods we calculate the 
asymptotic behaviour of the probability that no pairs of walkers from different families have met up 
to time t. For d > 2, this is constant, but for d < 2 it decays as a power t~", which we compute to 
0{e^) in an expansion in e = 2 — d. The second order term depends on the ratios of the diffusivities 
of the different families. In two dimensions, we find a logarithmic decay {\nt)~°', and compute a 
exactly. 

1 Introduction 

Consider the following problem: N random walkers set off from the vicinity of the origin, in d- 
dimensional euclidean space, at time t = 0. They are divided into p different families: the number 
of walkers in the jth family is rij, so that N — X]j=i "-j- Within a particular family, walkers are in- 
different to each other: their paths may cross. However, each family behaves viciously towards all the 
others: if two walkers from different families meet, both are annihilated. We may ask many different 
questions about this problem, but a fundamental quantity is the probability P{{nj};t) that all the 
walkers have still survived up to time t. Equivalently, we may consider the ensemble of N indepen- 
dent random walks: P{{Tij};t) is the fraction of these in which no walkers of different families have 
intersected up to time t. 

For a discrete time process on a lattice, if rj^ (t) is the position at time t of the i^jth walker of the 
jth family, then P{{nj}; t) is the expected value of the indicator function 

n n n n (i-'5(r7(o,r-(0)) (i.i) 

t'=0 l<j<k<pyj = l Uk=l 

This problem is of interest for several reasons. It generalises a number of cases: 

1. rij = 1 (1 < j < p) corresponds to vicious walkers^ a term first introduced by Fisher |Q. It 
has been studied using e-expansion methods ||] similar to those of the present paper. The 
survival probability is known exactly for c? = 1 in the case when all walkers have the same 
diffusion constants: it decays as a power [|l], ^, ||. These methods rely essentially on a 
fermionic description of the problem ^ . Krattenthaler et al. [Q introduced the method of the 
Schur functions and Katori and Tanemura Q developed this and discussed the relation with the 
random matrix theory. These methods do not appear to extend to the case when the diffusion 
constants are different. Results in this case have been reported for p = 2 

2. The case p = 2, with ni = n and n2 = 1, has been studied by Krapivsky and Redner^, ||, ^ as a 
model of n predators ('lions') hunting a prey ('lamb'). They were able to obtain exact results for 
the asymptotic behaviour of the survival probability, again in d = 1, for the cases n = 1,2 and 
arbitrary diffusion constants. For general n, the exponent is related to the smallest eigenvalue of a 
Dirichlet problem in a certain (n — l)-dimensional compact region, and is probably not expressible 
analytically, but for large n these authors were able to estimate its behaviour. The 'lion-lamb' 
problem for d = 1 is related to a version of the 'ballot problem' in which it is required to know 
the probability that one candidate in a ballot remains ahead of the n others at all stages in the 
poll. Exact results are known only for n <2 [To[ |. 
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3. The 'lion-lamb' problem has another interpretation, in terms of multiscaling: if we first regard 
the trajectory £ of the lamb as fixed, and if Pi{t) is the probability that it has not been met by a 
single lion, then 

P{n,l;t) ^ {pe{tr)f (1.2) 

where the average is over all the realisations of £. The fact that P{n, l;t) decays with t with an 
exponent which is not simply linear in n is symptomatic of multiscaling in this problem. 

4. More generally, we can regard P{ni, n2, ■ ■ ■ , Up] t) as being the average of the nith power of the 
survival probability of a single walker of family 1, in the presence of (ri2, . . . ,np) walkers of the 
{p — 1) other families. 

5. Our problem has a strong resemblance to that of the intersection probabilities of Brownian paths. 
In this case, one studies the ensemble of N random walks in d dimensions each of which begin a 
distance 0{a) from the origin and which arrive on a hypersphere of radius r = R ^ a before they 
cross r = a, irrespective of how long this takes. Once again the walkers are divided into families, 
and in this case one is interested in the probability P{{nj},R^a) that the paths of walkers of 
different families do not intersect. Thus, instead of ( |1 . 1[ ) , P is the expected value of 

Uk oo oo 

n n n n n (i-'5(r;^(o,r^M^"))) (i-s) 

l<j<fe<p Uj=l Vk=l t'=Q t"=0 

and it is supposed to decay as (R/a)^" as R/a — > oo, where a depends nontrivially on the 
{rij}. This problem is trivial in d = 1, and turns out to have an upper critical dimension d — A, 
below which an e-expansion is possible jl^. For d — 2 a,n exact formula for a{{nj}) has been 
derived |l^, by exploiting the conformal invariance of the problem. 

Given these remarks, it seems important to investigate the general case described in the opening 
paragraph. As far as we know, the fermionic methods used to attack the vicious walker problem for 
d = 1 do not extend to this case. We have therefore employed a renormalisation group (RG) method, 
which yields, for d < 2, results for the exponent a{{nj}) of the power law decay of P{{nj};t) as a 
power series in e = 2 — d. By using field-theoretic methods, the calculation is streamlined, and, once 
the formalism is set up, involves relatively little explicit calculation. We have carried this computation 
through O(e^), and for arbitrary diffusion constants of each family. It would be tedious, but not difhcult, 
to carry it further, as the actual Feynman integrals are elementary. We also show that in two dimensions 
P{{nj}; t) decays as a universal power of Int. 

The layout of this paper is as follows: in Sec. ||, for completeness, we collect all our results and show 
how they reduce in the above-mentioned special cases. In Sec. ^ we set up the field-theoretic formulation 
of the problem, then in the next section carry out the RG analysis. Sec. |] contains a summary and 
further remarks. Several of the detailed calculations are relegated to Appendices. 



2 Results 



Let p be the number of families, rij be the number of walkers in the jth family, and Dj be their 
diffusivity. Let P{{nj};t) be the survival probability 



P{{n,};t)^E 



n n n n {i-Kr7in<^m 

t'=0 l<j<k<p Uj = l i/fc = l 



(2.1) 



2.1 d>2 

In this case there is a finite probability that any pair of walkers will never meet. As a result, P({nj}; t) 
approaches a non-universal constant value less than 1, with leading power-law corrections of the form 

l(2-d)/2_ 
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2.2 d<2 

P({nj};<) - const. as t ^ oo, (2.2) 

where, with e — 2 — d, 

a = J^is + Ti £^ + Oif") (2.3) 

with 

1<J1<J2<P 

= 3(C?-C2), 

n^j rij^ rij^ 

l<il<j2<j3<P 

+ 3 E { ('^ii - 1 ) ^(^ii ' > ^ii ) + ("i2 - 1 ) In ^(^i2 . ^ji > ) } 

1<J1<J2<P 

'i<ik</' V (i^..+i^.2)^ j+".2ln(^ (i^,,+Z),J^ 

- 1 V n n i^/ ^ji-Pi2p.-.+2g,J(2jJ,,+I),J \ 

l<3l<h<P ^ \ 01^ 32) / 

where 

Cfe = Es' fc=l,2,..., (2.5) 

and 

Prom this may be deduced various special cases: 

2.2.1 Equal diffusion constants 

Assume that Dj = £) for all j = 1,2, - ■ ■ ,p. Then 

a = J^ie + J^2ln^s^ + 0{e^) (2.7) 

with 

i<ii<i2<p 

i<ii<i2<i3<p i<ii<i2<p i<ji<i2<p 

= H^i - ^1 - 2^?i<^2 + C2 + C3) , (2.8) 

Note that these are expressed in terms of symmetric polynomials in the {rij}. This in fact holds to 
all orders in e. 
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2.2.2 Vicious walkers with unequal diffusion constants 

When 

I otherwise, (^■^) 

a should be equal to the survival exponent ijjs.p of the vicious walkers. In this case Ck = p for k — 1, 2, • ■ • 
and the result (|2.3D gives 



'j=l (1<J<P), nfc=0 (fe>p+l) 
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(2.10) 



2.2.3 Vicious walkers with equal diffusion constants 



The result gives 

V'S,p = 0:\nj = l (l<j<p}, nk=0 (fe>p+l) 

= ^(2)^ + 1(3) In I + 

= \p{p~l)e+\p{p-l)ip-2)\n^e^ + Oie'). (2.11) 

This agrees with the result reported as Eqn. (5.2), with Eqn. (3.13) in Mukherji and Bhattacharjee|Q 
(see also [||). 

It has been proved that [0, ^, ^ 

V'S.p = iP(p - 1) for d=l {i.e. e=l). (2.12) 



Note that although this exact result agrees with that from the first-order e-expansion (2.11) on setting 
e = 1, this is probably fortuitous, as, in the case of unequal diffusivities the exact result depends on 
their ratio, while the first-order term in (2.1C) does not. 

2.2.4 'Lion-lamb' problem with unequal diffusion constants 

The 'n lions and one lamb' problem studied by Krapivsky and Rcdner |Q, ^ is a special case of the 
present model in which 

n for J = 1 

1 for j = 2 (2.13) 
otherwise. 

In this case Ck — I + n'^ for k — 1, 2, ■ • ■ and the result ( pT^ ) gives 

~ Q^lni=n, 712 = 1, ^3=0 {j>3) 

= ine+in(n-l)ln(^^i±^) e' + 0{e'), (2.14) 
where rj = D2/D1. Redner and Krapivsky |^ reported the exact solution for n = 2 in d = 1 (i.e., 

'n 



IT I + rj 

It was shown that p^^^"^^ (rf) is monotonically decreasing in 77 and 



(2.15) 



/3r'*(0) = 1, I3r^\l) - I, hm /3r'*(?7) = \. (2.16) 
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If we neglect 0{e^) and set n = 2, e = 1 in (2.14), we have 



/3— (^) = l + iln(^i 



l + 2?7 



-77)2 



which is monotonically decreasing in rj and 



„ox.^Q^ = 1, i3^PP'°-{l) = 1 + i In I ~ 0.856, lim fS^^'^'"''- (r,) = -oo. 
2.3 Two dimensions 

In this case, there is a logarithmic decay with universal exponent: 

- const. (In t)-" + O (^-^^^ (2.17) 

where 

1<J1<J2<P 

Note that this is independent of the Dj (as long as no pair of them both vanish) : the dependence shows 
up only in the prefactor and the non-leading terms. 



3 Field-theoretic formulation 

In this section, we set up the general problem as a continuum field theory, so that the powerful techniques 
of the field-theoretic RG may be applied. The general method for formulating such stochastic particle 



systems as field theories, as originally proposed by Doi|14| and Peliti||l5|, has been described at length 
elsewhere and we shall only summarise how this is applied in the case of interest. 

Initially, the problem is formulated on a lattice, for example Z'', the sites of which are labelled by 
a vector r. The microstate of the system at a given time is specified by occupation numbers {TOj(r)}, 
which specify that there are mj{r) walkers of family j at the site r. Note that we treat walkers of 
the same family as identical particles: this makes no difference in the problem of interest. To each 
microstate is associated a vector in a Fock space J-, built by applying raising operators to the vacuum, 
or empty state, |0): 

i{-.(r)})= n n«](rr^""'io> , (3.1) 

where [aj(r), aj^(r')] = SjkSrr', and aj(r)|0) = 0. Let p{{mj{r)}]t) be the probability of finding the 
system in this microstate at time t, and define the state G J-^ 

!*(<)) ^ p{{m,ir)};t)\{m,{r)}) . (3.2) 

{mj{r)} 

Then the master equation, which is linear equation describing the time-evolution of the probabilities 
p({mj{r)};t), is equivalent to the Schrodingcr-like equation 

d|*(<))/fit = -i?|1'(t)) (3.3) 

where H : J- ^ J- may be expressed explicitly in terms of the raising and lowering operators. For the 
case of independent random walks in continuous time, 

H^H,^j2{D,/b') J2 (a](r)-a](r'))(a,(r)-a,(r')) (3.4) 

J = l (r,r') 
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where b is the lattice spacing, and the sum is over nearest neighbour pairs of sites (r,r')- 

The probabihty of finding the walkers at sites rj^ at time t (where 1 < i^j < nj with 1 < j < p) is 
then given by 



l*(0)) 



(3.5) 



j=i vj-- 



Of course, when this is summed over all the rj^ , it gives unity. 

Before considering how to implement the non-intersection constraint, let us first discuss the con- 
tinuum limit and the path integral representation. In this non-interacting case, the continuum limit 
may be taken rigorously. The raising and lowering operators go over into (distribution-valued) field 
operators satisfying [(/)j(r), (j)\{v')] = 5jk6{r — r'), and the generator of time evolution becomes 



(3.6) 



Since the walkers are all supposed to begin in the vicinity of the origin at t = 0, that is, a finite number 
of lattice spacings away, in the continuum limit 6 — s- 



|vl/(0)) = 0t|0) EE f[(0t(o))"^ |O) 



(3.7) 



The path integral representation is derived||l6| by breaking the time interval (0, t) into slices of 
length At, so that the time-evolution operator e~*^" is the product of factors e~^*^" « 1 — AtHo, 
and inserting a complete set of coherent states at each time slice. This has the effect of replacing the 
operators (j)j (r) and 0j (r) by time-dependent c-number fields {t, r) and (j)* {t, r) respectively. After 
taking the limit At — > 0, the matrix element (|3.5|) becomes a functional integral 



where 



and O* = 11^=1 ' 



j = l j = lUj=l 



So 



dtd'^r 



(3.8) 



(3.9) 



3.1 Interactions 

We now discuss how to incorporate the constraint that walkers of different families should not meet. 



Rather than insert the indicator function (1.1) into the path integral, it more convenient to consider 
a slightly more general problem in which, before taking the limits 6 — )■ and At 
trajectories is weighted by a factor 



0, each set of 



nn-P 




JU2 



(6VAt)m,,(t',r)m,,(t',r) 



(3.10) 



where the Xj-^j^ (b'^/At) > are a set of dimensionless parameters (the factors of b and At are inserted 
to make the continuum limit simpler.) The case of strict non- intersection corresponds to the limit 
Xjij2 — > oo. However, we shall show that, for d < 2, the leading behaviour is independent of the precise 
value of these parameters (as long as they are all strictly positive) and, moreover, the RG fixed point, 
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(s + Djk^J^ 




-X 



hh 



Figure 1: Feynman rules for the interacting theory ( 3.11 ). Time always flows towards the left 



at which non-leading corrections to the asymptotic behaviour disappear, corresponds to the limit of 
infinite Xj^j2 ■ 

In the formal continuum limit, this corresponds to a modification of the action in the path integral 

A,,,. / ^*6l(l)j,(l)j,dtd'^ (3.11) 



S — So + ^ ^Vjij2 

i<ii<i2<p 



3.2 Feynman rules 

The Feynman rules for this theory are very simple[^ and are illustrated in Fig. [|. We denote averages 
and correlations with respect to the bare action 5*0 by the subscript o- averages with respect to the full 
action S are denoted by (•). 



• the Fourier-Laplace transform of the bare propagator 



Gf'')(s,k)o 



dt I d'^re-'^e'^ " 



(0,(t,r)0*(O,O))o 



(3.12) 



is represented by a line directed towards increasing time (conventionally, right-to- left). In the 
(<, k) representation the bare propagator is simply e~^^^ *; 

• the interaction (— Aj^j^) is represented by a vertex with one incoming and one outgoing pair of 
lines each of type ji and j2 ; 

• As usual, wave number k and (imaginary) frequency s are conserved at the vertices, and internal 
loop integrations J{ds/2TTi) and J {d'^k / {2tt)'^) are carried out. 

3.3 Renormalisation and operator product expansion 

The survival probability is now given by the correlation function 



Go{t) 



j = l Vj=l 



(3.13) 



evaluated with the weight e"'^. 

However, this does not exist in the formal continuum limit, because the perturbative Feynman 
diagram expansion of G contains ultraviolet (short-distance or short-time) divergences. Physically 
this is because two walkers, having interacted once, are, in the continuum limit, likely to interact 
an infinite number of times as At 0. This divergence may be regulated, either by imposing an 
explicit cut-off < A in the Feynman integrals, or, more easily, by dimensional regulation. For 
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d < 2 this field theory is renormalisable: the singular dependence on the regulator may be absorbed 
into a finite number of parameters. In the case of the theory of interest, this procedure is particularly 



simple! 16 : no renormalisation of the field (f>(t, r) nor of the diffusion constants is required, only a simple 
renormalisation of the coupling constants Xjji , which can be computed exactly to all orders. The lack of 
field and diffusion constant renormalisation holds mathematically because there are no loop corrections 
to the propagators. Physically it is because an isolated walker does not interact, even with itself, in the 
absence of any branching processes. When the coupling constant renormalisation is done, all correlation 
functions of products of ((> and 0*, at distinct space-time points, have a finite limit as the regulator is 
removed, when expressed in terms of the renormalised couplings. The fact that the renormalised theory 
must be defined at some arbitrary scale then leads to RG equations. 

However, this procedure is not sufficient to render finite correlation functions involving so-called 
composite operators like O* = YVj^i ^ji^'^)"^' ■ Physically, this is because if the walkers all begin 
at exactly the same point, they will all annihilate each other immediately! In order to obtain finite 
renormalised correlation functions, it is first necessary to point-split the fields: 

J]^*(0,or -^n0;(O,r,r^- (3.14) 

(Note that it is not necessary to split the starting points of walkers of the same family, since they do 
not interact.) Now consider a correlation function of this product with an arbitrary product A of fields 
whose time arguments are all strictly positive: 

{Af[cb*iO,r,D (3.15) 

In the cut-off theory, we could simply make a Taylor expansion of this in powers of the r j . This would 
have the form 

{AO*{0, 0)) + Cn{{r,}){AO:{0, 0)) (3.16) 

n 

where the summation is taken over all possible derivatives {O,*} of Y[j 4>j{0,J^j)"'^ with respect to the 
{rj}. On the basis of dimensional analysis, the first term gives the lea ding behaviour for Ga{t) as 



t — > cxD, at least at c? = 2. In the interacting theory, however, each term in ( 3.16 ) has to be renormalised 
separately. As a result 

{Af[cb*{0,r,r^) = Z-\AO*j,)+Y,ZalCn,Ri{r,}){AOlR) (3-17) 

j = l n 

where all the correlation functions are finite as the regulator is removed. Since each of the {On,B.} 
may acquire a nontrivial anomalous dimension through this procedure, the renormalised functions 
C,i.fl({rj}) have a non-trivial dependence on their arguments. The important feature of ( 3.17| ) is that 



the renormalisation constants Zq, Zq^ are independent of A. This we may write the operator product 
expansion (OPE) 

n</.*(0,r,)"^ = Z^'OU0,0) + J2Zala,M{r^,})O:^j,i0,0) (3.18) 

j=l n 

where the OPE functions Crt.j;({r,-}), etc. are in general nontrivial. 



Each term in the OPE ( 3.18 ), when substituted into G, will give rise to nontrivial power-law depen- 
dence on t for d < 2. However, since in the noninteracting theory we know that only the first term is 
important as t — > oo, we shall assume that this remains true for sufficiently small e. Further discussion 
of this point will be postponed to Sec. ||. 
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Figure 2: Some of the diagrams contributing to the survival probabihty Go- The case ni = 2, 712 = 1, 
713 = 1 is shown as an illustration. The one-loop diagram shows an interaction between a walker of 
family 1 with one of family 2, and is proportional to A12. 



3.4 d>2 

In the absence of any interactions, the survival probability Go{t) = 1, as can be seen by evaluating the 
first diagram in the expansion shown in Fig. || with all the wave numbers qj^ set to 0. For d > 2, the 
higher order terms give contributions which correct this constant, as well as terms which are subleading 
as i — > 00. For example, a typical one-loop diagram like that in Fig. ^ is proportional to 

- J dt' j d'^fce-(^-+^-)'' * oc - y ^ ^ '- d-'k (3.19) 

This integral diverges at large |k| for d > 2, a consequence of taking the naive continuum limit. If we 
impose a cut-off |k| < A = 0(6), the leading term behaves as K'^~^ , with a non-universal coefficient, 
while the remainder is finite as A ^ 00 and behaves as t"^'*"^)/^. The non-universal constant term 
corresponds to a finite probability that the walkers annihilate at short times, before escaping each other. 
This behaviour persists at higher orders in the interactions. 

For d < 2, however, each successive term in the bare perturbation grows as a larger and larger 
positive power of and it is necessary to resum the expansion. The renormalisation group (RG) 
provides a consistent framework within which to carry this out. 



4 Renormalisation Group Analysis 
4.1 Coupling constant renormalisation 

As usual, the renormalised couplings Xrju are defined as the values at the normalisation point of the 
irreducible vertex functions F^^'^^ ((s^-, q^), (s^,q^); (sj,qj), (sk,^k)), which are the truncated Laplace- 
Fourier transforms of {(j)j{t'j,rj)(j)j{t'i., r'i^)(j)*{tj,rj)(j)*{tk, r^)). (There is no field renormahsation in this 
theory.) It is convenient to choose the normalisation point as 

4 = s'k = = Sk=(T ^0; q^ = q^,. = q^ = q^ = (4.1) 

This class of theories has the special property that the renormalised coupling constants may be computed 
to all orders JT6|. The calculation is summarised in Appendix A. The result is 

where 

- ^5|^r(l - dl2) = l/{4n) + 0{e). (4.3) 
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It is convenient to define the dimensionless renormalised couplings as 

2 \ (i/s 



gR,k^\R,k (2'^)"'^' ■ (4-4) 

As will be seen, these are the natural expansion parameters for the renormalised perturbation expansion. 
4.2 Renormalisation of O* 

As discussed in the previous section, we are interested in computing the asymptotic behaviour at large 
t of 

GoW = / n n rf'V?(0.(i,r7)O*(O,O)) , (4.5) 

in the regularised bare theory, where O* = 11^=1 i'Pj)"' ■ However, for the purposes of renormalising 
O* , it is more convenient to choose the time arguments of the fields to be independent, and to 
consider the laplace transform with respect to these times. Thus we define 

Goiit"/}) = / n n dV/(q^,{f/y/)0*iO,0)) , (4.6) 

j=lu,=l 

and 

= / n n } Go({t?}) (4.7) 

From this we may define the irreducible vertex function [pT| , by truncating the external propagators: 

ra({.r};{A,4) = (4.8) 

n^.rn::;.r(^r) 

The renormalised vertex function is 

ToR{{s''/h{gR,k},<j) = Zo({A,fe},a)ro({s7};{A,4) (4.9) 
where Zq is fixed by the normalisation condition 

roi?({s;n; W},^)Li^ = 1- (4.10) 



In writing ( |4.9| ), we have made it clear that that the (un)renormalised vertex function is to be thought 
of as depending on the (un)renormalised couplings. 

Now, although we have used the condition ( |4.1C| ) on Tqh to define Zq, the same multiplicative 
renormalisation also renders 

GoRit;{XRjk},<j) = ZoGoit;{Xjk}) (4.11) 



finite, for t > 0, where Go is defined in (4.5). For this to be true, it is important that the fields (pj {t, r'^' ) 
are not evaluated at the same point. This would lead to further UV divergences. However, these occur 
on a set of measure zero in the integration in (4.5), and are harmless. 

4.3 Callan-Symanzik equation for Go 

Define the RG functions 

P,k{9Rjk) = a(^^) . (4.12) 
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ro = 




Figure 3: One-loop diagrams for F 



and 



lo{{9R]k}) = ( (J—\nZo{{\jk},(T) 



The fact that a{d / d<j)Go\{\jk) — then iniphes the Callan-Symanzik equation 



(4.13) 



d \ — 

GoR{t,{gRjk},cr) = 0. 



l<j<k<p 



dgRjk 



(4.14) 



If the couphngs {gRjk} flow towards a nontrivial fixed point {<?|jjj,} at which Pjkigjijk) — i^^ '^6 
shall show happens for d < 2), then in estimating the leading asymptotic behaviour as ct — > oo it is 
sufficient to replace (114) by 

d 



da 



ih GoRit,(T) = 0, 



(4.15) 



where = JoiigRjk})- This has the solution GoRit, cr) oc cr^o, as cr ^ cxd at fixed t. 

However, simple dimensional analysis implies that GoRit,(y) is function of only the combination 
{(jt). Hence we find that 

GoRit, cr) ~ const, t"" (4.16) 



with 



a 



-Jo 



(4.17) 



4.4 /3-functions 



From (4.2) and (4.4), we find after some algebra (see App. A) 

PjkigRjk) = -\{£gRjk - bdgjijk). 



(4.18) 



Note that this is exact to all orders in gRjk, and that, for e > 0, there is an infrared stable fixed point 
at 

g*H,k^e/b,i = 47:e + Oie^) , (4.19) 
whose value is independent of the diffusion constants. 



4.5 One-loop calculation of 70 



Consider the expansion of the vertex function (4.8) as a power series in the coupling constants {Xjk}- 
To first order, this is given by the sum of diagrams like that in Fig. |3| explicitly 



To = 1 - ^1 ('^ii ' ; ^ii ' -Dj2 ) + ( { A^j^ } ) , 

1<J1<J2<P 



(4.20) 



where the integral Ii{Sj_^, Sj^; Dj-^, Dj^) is the same as occurs in the coupling constant renormalisation: 
see Appendix A. The combinatorial factor nj^rij^ counts the number of ways the walkers from family 
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ji can interact just once with those of family j2. The renormahsation constant Zq is then the inverse 
of this evaluated at the normalisation point s'^^ — a. Thus 



\nZo= ^ 

l<jl<j2<P 



d/2 



(4.21) 



and so, by (4.13) 



70 



32 



d/2 



i<jl<j2<P 

'\ XI %l"j2 X Mi?,JlJ2 +C'({5flJfc})■ 

l<jl<i2<P 



Next we set gnj^.j^ = g*. By (4.1£) 



i<ii<j2<p 



(4.22) 



(4.23) 



Through (4.17), this gives the result (2.3) up to 0{e). We remark that to this order the result is 
independent of the diffusion constants, as long as no pair Dj^ + Dj^ vanishes. This last is of course a 
pathological case, since then the two families are immobile and cannot meet. 

4.6 Two-loop calculation 

There are three types of diagrams contributing to To at order A^. They are illustrated in Figs. ^-c. 
Contributions of types (a) and (b) involve the one-loop integrals Ii defined above, while those of type 
(c) involve 



{(■'^^2 + 4) + (Dj2 + ^i3)(k + q)2}{(5;, + sj, + 4) + Dj^k^ + Dj^cf + D,,{k + q)^} ' 



(4.24) 



Define 



(4.25) 



Then 



^ — 1 - ^ Hj^ Uj^ Xj^j2 h (cr; Dj^ , Dj^ ) 

l<jl<j2<P 

+ '^ii"i2(Aiu2)^(A(cr;Dji,LijJ)^ 

l<jl<j2<P 



+ %1'^J2"J3"-J4 {^jU2-^j3j4A(CT;£'ji,£'j2)A(CT;L'j3,L'j4) 

l<jl<j2<j3<j4<P 

j3 '^j2j4 A (o'; ; ^j3 ) A (c^; ^j2 1 ^j4 ) + -^ji j4 •^j2j3 A (ct; i^ji , Dj^ ) A (ct; i'ja : ^j3 ) } 

+ H {"jlKl - 1)"-J2%3^J1J2^J1J3A(ct;£'ji,£'j2)A(ct;£'ji,£'j3) 

l<il<j2<j3<P 
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(a) 




+nn (jij^ - l)nj3 A^.j^ h {a; Dj^ , Dj^ )/i {a: Dj^ , Dj.^ ) 
i ^ rij, (uj, - l)nj., {rij., - 1) (/i (cr; Dj^ , D^Jf 

1<J1<J2<P 



+ X! '^■Ji '^^3 { '^Ji J2 j3 ^2 (o-; IJji , , -Dj3 ) + Aj, j2 Aj, j3 /2 (cr; ZJj^ ; . -Djs ) 

i<ji<i2<i3<p 

+-^jij3-^j2j3-^2(cr; Dj-^,Dj^,Dj^) + ^^jijaXj^j^hicr; Dj^,Dj^,Dj^) 

+ ^3233 ^3133 ^2 (cr; £'j2 i ^^33 ' -^jl ) + ^j233 ^3l32 h {<7\ Dj^ , Dj^ , L»j ^ ) j 

+ XI '^Ji ('^ii ~ (•^jii2 )^-^2 (cr; Dj, , Dj^ , D^, ) 

l<jl<j2<P 

+ E %i^J2K2 - l){\nj2fh{'y: D,,,D,„D,,) + Oi{X%}) (4.26) 

l<jl<j2<P 

Each term in this sum corresponds to a diagram of class (a) to (c3) in Fig. |^. The combinatorial factors, 
polynomials in the rij, count the number of ways different walkers from a given family can contribute 
to each of these processes. (It is simplest to check these factors for small values of p and iij.) From 
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Appendix A 



h / 2 \ ''^^ 

A(.;I.„i..) = f U-^ (2.r/^ (4.27) 



Moreover, as shown in Appendix B 



(i/2 / r, \ d/2 



Dj+DkJ \Dk+D 



l^-l-lnR{Dj,Dk,D,) + 0{l) 



bl{2a)-', (4.28) 



where bd is given by (4.3) and R{Dj, Dk, Di) is given by (2.6). 

Next we compute In through 0{{X^f.}), perform the differentiation ad /da at fixed bare coupHngs 
Ajfe, then re-express the resuh as a series in the gRjk, to the same order. The result has the form 

lo{{9Rjk}) = Si + Ji?2&^ + i?3 + 0({<?|, J) (4.29) 

with 

i<ii<j2<p 

l<il<j2<P 1<J1<J2<P 

+ X 'nn^h'^33^34i9Rjij29R33U + 9Rjij39R32H + 9Rjiji9R]2]3) 

+ X TT-jinnnj,[inj^ - ^)9Rnj29Rjij3 + ("is ^ ^)9Rjij29Rj233 

I<jl<j2<j3<p 



+ i'^j3 - ^)9Rjlj39Rj2j3} 



+ ^ X "jlKl - 1)"J2(«J2 - l)(5fljU2)^ 

i<ji<i2<p 

+ X "ji"-j2'^j3(5i?jii2 5flj2j3 + 9Rj2j39Rji33 + 9 R jijs 9 R 3132) 

1<J1<J2<J3<P 

+ 5 X '^Ji"j2('^ji +"j2 -2)(ffi?jij2)^ - n X ^3i^h9Rjij2\ > (4-30) 
i<ii<j2<p \i<ii<i2<p / 

X "jl%2"j3{5KjU25flj2j3lni?(£'jl,-Dj2:-Dj3) 
l<il<j2<j3<P 

+9Rj2j39Rjij3 ln^(^j2i-Dj3,L'j J + 9Rjij39Rnn Ini?(-Dj3, -Dji, -Dja)}^^ 

-3 X '^Ji'^J2{("ji - l)lni?(L'ji,£'j2,i:'ji) + (rij, - 1) lni?(L>j, , Llj^ , j|fe^(.gflj,j 
i<ii<j2<p 

An important check of this calculation is that the double poles in e in the two-loop contributions 
are cancelled by the coupling constant renormalisation in the one- loop contributions. This has the 
consequence that ^o{{9Rjk\) has a finite limit as e — » 0, that is, B2 vanishes. This is shown in 
Appendix C. 

Now we set 

9R'jk =g*R = — for aU (j, fc), (4.31) 

Od 



following (4.19). Then we have 

Yo = Ble + Ble^ + 0{e^) (4.32) 
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with 



' 2 5Z "'■'1 

i<ii<j2<p 



{ In R{Dj^ , Dj^ , ) + In R{Dj^ , Dj^ , Dj, ) + In R{Dj.^ , Dj, , D^, ) | 

l<il<j2<j3<P 

3 J2 { (%i - 1) In RiDj, , , Dj, ) + (n,, - 1) In , , D^, ) } 



i<ji<i2<p 



l<il<j2<j3<P 



"ii"j2"j3ln 



1<J1<J2<P 



3 E "Ji"j2ln 

l<jl<j2<P 



Through (4.17), this gives the resuh (2.3). 



4.7 Two Dimensions 



When d — 2 the couphngs {gn^k} ^I'e marginally irrelevant, that is they flow logarithmically slowly 
towards zero. In that case, it is not sufficient to set them equal to their fixed-point values, but instead 
the full Callan-Symanzik equation ( |4.14| ) must be solved. Using the fact that GoR{t, {dRjk}, f ) depends 
on (7 and t only through combination (cri), this may be rewritten 



d 



d \ ^ 



^di ~ ^'^^^9Rjk}) + f^jkigRjk) _ ) GoRit, {gRjk}, cr) = 



dgRjk 



(4.33) 



l<j<k<p 

The solution by the method of characteristics is standard. Define running couplings {gjk{u)} by 

'^^9jk{u) = -I3jk{gjk{u)) 
with initial conditions 5jfc(l) = gRjk- Then 

GoR{t, {gR,k}, a) = e/r 7o(fe..})(<i«/«) GoR{a-\ {g^'^t)}, a) 



In our case, 



9jk{u) = 



1 



\b2\nu + gj^^jf^ b2lnu 
so that, using (4.22) the exponent in ( 4.35D is 



+ C'((lnu)-2) 



InHat) J2 "Ji"j2+0((ln(at))~i) 
i<ji<i2<p 



(4.34) 

(4.35) 
(4.36) 

(4.37) 



Exponentiating this yields the result quoted in ( 2.17 ). The dependence of the last factor in ( [4.35 ) on 
{gjki'^t)} also generates corrections which are down by 0{{\n{at))~^ . All the non-universal behaviour 
resides in these, and higher order, corrections. Note the absence of c orrec tions 0(lnlnt/lnt), which 
may be traced to the lack of higher-order terms in the beta- functions ( 4.18 ) . 
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5 Discussion 



We have presented a generalisation of the vicious walker problem in which walkers from different families 
annihilate on meeting, but walkers from the same family ignore each other. We have studied the problem 
in a field-theoretic renormalisation group framework, suitable for understanding universal quantities 
such as critical exponents. We have focussed on the probability that all walkers have survived up to 
time t, and we have showed that, in dimension c? < 2, this decays as a power, t~"(f"3l), where rij is 
the number of walkers in the jth family. While this result is true to all orders in e = 2 — d, the actual 
values of the exponents can, by this method, be evaluated only as a power series in e, which we have 
carried out to second order. The coefficient of the 0(£^) term depends on the ratios of the diffusivities 
of each family, as well as the {nj}. The lack of dependence of the 0{e) term on this ratio may be 
traced mathematically to the fact that the same bubble diagram enters into the coupling constant 
renormalisation (Fig. H) and the renormalisation of the composite operator O* (Fig. ||). 

For the same reason, the exponent for = 2 walkers does not depend on the ratio of their diffu- 
sivities, because in this case the one-loop result is correct to all orders (as long as we do not expand 
the coefhcients in powers of e.) The same would be true, to first order for N > k, and to all orders for 
N = k, ii the two-body interactions we consider were generalised to fc-body interactions with fc > 2, 
although in this case the upper critical dimension would be reduced to dc{k) = 2/{k — 1). 

The value of a{{nj}) was shown to be related to the anomalous dimension of a certain composite 
operator C{„^j — Hj <^j"^- This structure means that our results can be straightforwardly extended to 
other physical observables. For example, the reunion probability R(t) that all the walkers have survived 
up to time t and are all located within a distance 0{b) ^ t^^^ of each other (where b is for example the 
lattice spacing), is related to a correlation function 

i?(t) - /(O(t,r)O*(0,0))d'V (5.1) 



where O = Y[j 4>]' ■ Since the theory is symmetric under {t —t, (f)j 0*), O has the same anomalous 
dimension as O* . Hence 

R{t) ~ const. i-(JV-i)(d/2)-2«({„,}) (5 2) 

where the first term in the exponent comes from simple power counting, and the factor of 2 in the 
second reflects the important fact that the anomalous scaling of composite operators at different times 
(and points) is multiplicative. 

The fact that C{„ .j is symmetric under permutations of the families j = 1, . . . ,p has the important 
consequence that the exponents a{{nj}) are also symmetric. The form of the e-expansion implies that 



this is true to all orders in e. It may be traced to the operator product expansion ( 3.1g ): the next-to- 
leading terms on the right-hand side, which do not vanish on integration over the spatial coordinates 
must contain at least two derivatives, for example 

(V0i.v02)0r"V'r"' n 'f'T (5-3) 

3<j<P 

Power-counting shows that, at d = 2, the contribution of such terms is at least 0{t^^) down on the 
leading term, and therefore, for sufficiently small e, they yi eld o nly corrections to the leading behaviour 



which we have computed. However, since each term in ( 3.18 ) is renormalised separately, each gives 
rise to an independent scaling exponent. It is a very interesting question whether the first term in the 
OPE is dominant all the way down to d = I. Such a result would imply that the asymptotic exponents 
(but not necessarily the prefactors) for the cases {ni, n2, na} = {2, 1, 1} and {1, 2, 1}, for example, are 
equal. Yet, in one dimension, the two problems are certainly not isomorphic, because the ordering of 
the families along the real line makes a difference. A very similar situation occurs in the problem of 
intersections of families of Brownian paths in two dimensions: in that case the exponents have been 
computed exactly[^2[ |l3| and are known to be symmetric, despite the fact that the ordering of the 
families around the annulus is a priori relevant. In this example, this symmetry is also suggested by 
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the e-expansion below the upper critical dimension, in this case dc — 4, pl| . However, a recent result 
of Bray and Blythe|lJ] suggests that the situation in d = 1 is not so straightforward for our problem. 
In their Eq. (4) they report a result for the survival exponent for a single lamb, with diffusion constant 
D' , with Nl lions to the left and lions to the right. The diffusion constant of the lions is D, and 
their result is an expansion to first order in D' / D. It is reported to depend on the asymmetry — 
as well as the sum Nl + Nr, while our analysis shows that it should depend only on the sum, to all 
orders in the e-expansion. 

There are at least two possible resolutions: 

1. the qualitative conclusions of the e-expansion break down in d = 1, cither because the non- leading 
terms in the OPE dominate, or through some more systemic failure; 

2. for d = 1 there is a qualitative difference between the case when walkers if different families 
strictly cannot pass each other, and that, more appropriate to the field theory approach, when 
the annihilation rate is finite and therefore the order along the real line is not preserved. 

In any case, it may be shown, by generalising the arguments of Krapivsky and Redner[^ ||, |^, that 
for infinite annihilation rate in d = 1 the exponents (and indeed all the non-leading exponents), 

are simply related to the eigenvalues of the Dirichlet problem in a certain compact region consisting 
of a spherical hyperpolygon on the sphere Sn~2- Thus, if the postulated symmetry were to hold, it 
would imply that the Laplacians in different regions which are related by permutations of the {fij} are 
isospectral. We plan to address this question in a future publication. 
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A Coupling constant renormalisation 

The diagrams contributing to T^f^'^J are shown in Fig. ||. In (s, q) space they give a geometric sum. We 
have therefore the exact result 

((4. , J , (4 ' J q.. )) = t-t — r .^/ ''^'^ i — ( a. i) 



where s = s'^^ + s'^^ = sj^ + Sj^ and q = q^^ + q'j^ = qj^ + qj^ . The renormalised coupling is the value 
of this at — = sj^ = sj^ = cr and q^-^ = q^^ = q^^ = q^^ = 0. Thus 

^Bnn = ^ (A.2) 

1 + XjiiJi{'^'-:Dj^,Dj^) 
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where 



Ii{a;Dj,,Dj.,, 



{2nY 



r(i -d/2) f 2 



d/2 



d/2 



The dimensionless couphng (4.4) is then given as 



9Rjij2 



d/2 



{2a) 



(2a) 



-e/2 



-e/2 



d/2 



(2a) 



-e/2 



Thus 



^3132 



-d/2 



1 - 9 R 3132 — 



(A.3) 
(A.4) 
(A.5) 



(A.6) 



(A.7) 



Differentiating this equation with respect to a at fixed {Xj-^j^, Dj-^, Dj^) and using the definition of the 
beta- function ( 4.12| ) 

i-e/2)^^^Mll^ = (- ^—^] (3,,,,{gB.n32) (A.8) 



1 - 9Rji32 - 



^-9R3i32T 



That is, 



P3l32{9Rjij2) - -hi^9R3l32 - bd9Rjij2) 



(A.9) 



B The integral I2 

At the normahsation point, the two-loop diagram Fig. ^ leads to the integral 



h{<J-Di,D2,D^) 
dfiq f d'^k 



1 



{2Tr)d J {27T)d (2a + [D2 + D:i)\<.^) (3a + i:iiq2 + D2(k + q)2 + D^^-k^) 
As usual, the denominators may be combined using a Feynnian parameter integration over x: 

d'^q f d'^k 1 



Jq J (27r)^7 i^^r {{D2+D3)k^+x{Di+D2W + 2xD2k-(i) 
The wave-number integrals are now standard, and yield 

^2 - 7^^r(2 - d)iD, + D2)-'"\D2 + D^r'/^a'-^ x J 



where 



J= C dxx-'^'^{2 + xf-^ (l-- 



.XD2 



i?2)(i?2+i?3) 



-d/2 



(B.l) 



(B.2) 



(B.3) 



(B.4) 
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This has a simple pole at d — 2, arising from the end-point at a; = 0. However, we also have to extract 
the finite part. This may be done by writing J = Ji + J2 where 



Ji 
J2 



£ 
1 

dxx-'^/^ (2 + x) 





d-2 



d-2 



xDj 



The second integral is finite at d = 2: 



dx 
^ 



[Di + D2){D2 + D3 



xDl 



(Di + D2)iD2 + Ds) 
\nR{Di,D2,D3) + 0{e) 



where 



R{D^,D2,D^) = 



D1D2 + D2D3 + D1D3 



{Di + D2)iD2 + D3) 
Recalling the definition (^|^) of bd, we therefore find, after some algebra, 

<i/2 / r, \ d/2 



D1 + D2 



D2 + D3 



{2a)- 



1 1 

2£2 ~ 4e 



-d/2 



-)d-2 



0{e) 



\nR + 0{l) 



(B.5) 
(B.6) 



(B.7) 
(B.8) 

(B.9) 



which leads directly to ( 4.28| ). It should be remarked that a crucial simplification in this calculation 
arises because 

^^^"''^ -\{l + 0{e^)) . (B.IO) 



r(i-d/2) 



2 4 ' 



C Verification that 5o = 



By dH), 



B2-\ ^ n-j^nj^{gRj^j^f 

1<J1<J2<P 

Uj-^ Tlj^ Tlj^ Tlj^ 

{9Rjih9R3-i3i + 9Rjih9R323i + 9Rjij49R3233) 

l<il<j2<j3<j4<P 

+ X! ^3l^3'2'^33iP'h9Rjij29R3l33 +'n-329Rjij29R3233 +f^h9Rji'h9R3233 ) 

l<jl<j2<j3<P 

^ X! '^3l^32'^33{9Rjlj29Rjlh + 9Rjlj29Rj233 + 9Rjlh9R3233 ) 

1<J1<J2<J3<P 

+ 5 '^Jl"j2(5KjlJ2)^ - 5 "jl"j2Kl +"J2)(5/?JU2)^ + I Z ' 

i<ii<j2<p i<ii<j2<p i<ji<i2<p 

+ Z '^h'^j2'^33{9R3i329Rj2j3 + 9Rj2j39R3i33 + 9 R jijs 9 R 3132) 
i<ji<i2<i3<p 



'Jl "-^2 \yR3i32 



l<il<J2<P l<jl<J2<P 

E '^.■l".2n.3".4(5flj-,j25fl.3.4+5K 31 33 9 R 3234 + 9R3l3i9R32j3) 

'^<31<32<33<U<P 



2 I 5Z '^h''^329Rji32 

i<ii<i2<p 
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^<jl<j2<j3<P 

+^ J2 '^'ji^h(9Rjij2f - h i J2 ^hnj29Rjij2 

i<ii<i2<p \i<ii<i2<p 

0. (C.l) 
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